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Abstract. We investigate here the origin of the gas 
observed in extended emission line regions surrounding 
AGNs. We use the technique of calcium depletion as a 
test to prove or disprove the existence of dust in such a 
gas in order to discriminate between two main theories: 
(1) a cooling process from a hotter X-ray emitting phase 
surrounding the galaxy, (2) merging or tidal interaction 
between two or more components. We have obtained long 
slit spectroscopy of a sample of objects representative of 
different galaxy types although our main interest focus 
on radio galaxies. The spectral range was set to always 
include the [Call] AA7291, 7324 doublet. The faintness or 
absence of such lines is interpreted as due to the depletion 
of calcium onto the dust grains and, therefore, is a proof 
of the existence of dust mixed with the gas in the EELRs. 



1. Introduction 

Extended emission surrounding AGNs has been detected 
in many active galaxies. Such structures have spatial scales 
that vary from lkpc to hundreds of kpc. Many of the ex- 
tended emission line regions (EELRs) are found in ellip- 
tical galaxies which are not gas rich objects, thus: where 
does the gas come from? Does it come from stellar evolu- 
tion as a consequence of stellar winds and SN explosions? 
Or, on the contrary, is the origin external, and the gas 
is accreted from a companion galaxy rich in gas or from 
a hot halo surrounding the galaxy? In Paper I (Villar- 
Martm & Binette 1996) we made a review of the observa- 
tional evidences which supports or contradict these theo- 
ries, being an external origin easier to reconcile with the 
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observations. The morphological similarities with the sys- 
tems of filaments surrounding galaxies near the centers of 
rich clusters suggests that the origin of the EELR result 
from cooling of the hot gaseous halos which surrounds the 
central galaxies. On the other hand, morphological and 
kinematical studies of the gas often suggest the existence 
of a collision or merger in the recent past. 

The existence or not of dust mixed with the emitting 
gas has direct implications about the most plausible sce- 
nario for the formation of these extended ionized struc- 
tures. In the scheme of the cooling flow theory, the hard 
intracluster radiation field would prevent dust from form- 
ing. If the material consists of galactic debris, the dust/gas 
ratio is expected to have a value appropriate to the chem- 
ical composition of a normal galaxy. Our aim is to check 
if the dust does exist in EELRs. 

Ferland (1993) has proposed that the absence of the 
forbidden [Call] lines AA7291J324 (Fl and F2 hereafter) 
can be used to infer the presence of dust mixed with the 
gas in the Narrow Line Region (see also Kingdon, Ferland 
& Feibelman, 1995). Calcium is very sensitive to dust. It 
suffers strong depletion into the dust grains and there- 
fore, its abundance in the gaseous phase (responsible of 
the line emission) is lower than in a dust free nebula. The 
line emission should be fainter that expected, or even un- 
detectable. 

In Paper I, we proved that such test remains valid (and 
is very sensitive) under the conditions of the EELRs stud- 
ied here. We studied in detail all the most plausible al- 
ternative mechanisms to that of internal dust for explain- 
ing the absence of [Call] lines. No acceptable alternative 
solution was found and we concluded in favour of the va- 
lidity of the method initially proposed by Ferland (1993). 
The observational results and their interpretation are pre- 
sented in this paper. We have applied this test to a sample 
of objects of different types: radiogalaxies with extended 
emission line regions, cooling flows and starburst galaxies. 
Our main interest concerns to the EELR in radio galaxies. 
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Table 1. Observing Log 



Name 


RA(1950) 


Dec(1950) 


z 




Comment 


Exp time (s) 


PA 


NGC1052 


02 


38 


37.33 


-08 28 09 





.005 


Liner 


1800 


270° 


NGC6215 


16 


46 


47.0 


-58 54 30 





.005 


Liner 


2700 


240° 


NGC7552 


23 


13 


25.0 


-42 51 24 





.005 


Liner 


1200 


270° 


NGC7714 


23 


33 


40.59 


01 52 42 





.009 


Starburst 


2700 


270° 


PKS1404-267 (nuc) 


14 


04 


38 


-26 46 51 





.021 


RG 


2700 


270° 


PKS1404-267 (5" S) 
















2700 


270° 


PKS2014-55 


20 


14 


06 


-55 48 52 





.061 


RG 


3600 


190° 


PKS2152-69 (nuc) 


21 


52 


58 


-69 55 40 





.028 


RG 


2700 


270° 


PKS2152-69 (cloud) 
















3600 


290° 


PKS2158-380 


21 


58 


17 


-38 00 51 





.033 


RG 


2700 


270° 


PKS2356-61 


23 


56 


29 


-61 11 42 





.096 


RG 


3300 


285° 


PKS2300-18 


23 


00 


23 


-18 57 36 





.129 


RG 


2700 


240° 


2A 0335+096 


03 


35 


52 


09 48 10 





.035 


CF 


5400 


147° 


A2029 


15 


08 


30 


05 57 00 





.077 


CF 


2700 


270° 


A2597 


23 


22 


42 


-12 23 00 





.085 


CF+RG 


2700 


197° 



We describe the observations, data reduction and anal- 
ysis of the spectra in section 2. In section 3 we present 
the results of the comparison between measurements and 
model predictions. A detailed analysis of the nuclear spec- 
tra is presented in section 4 and conclusions compose sec- 
tion 5. 

2. Observations and data reduction 

The observations were carried out on the nights 21-23 Au- 
gust on 1993. All the spectra were obtained at La Silla 
Observatory, Chile, with the 3.6 m telescope, using the 
EFOSC 1 spectrograph with a CCD detector (TEK#26) 
of 512 x 512 pixels 2 of 27 ^m 2 . The slit width was 1.5". 
The grism used was, R150, with a dispersion of 120 A 
/mm, a wavelength bin of 3.3 A /pixel and covering a 
spectral range of ~ 6870-8560 A . The observing condi- 
tions where photometric. Table 1 gives a log of the obser- 
vations. 

2.1. Basic data reduction 

The reduction of the data was done using standard meth- 
ods provided in IRAF. The spectra were bias subtracted 
and divided by a flat-field frame (dome flat-field). Illumi- 
nation corrections along the slit were found to be negligi- 
ble. 

In general, we obtained three frames for each object 
and each slit position, allowing the direct removal of cos- 
mic ray events. For a given object, all the frames corre- 
sponding to the same slit position were averaged together. 
In the cases where only two frames were available, the 
cosmic rays were removed visually, replacing the affected 
pixels by the mean of the surrounding region. The spectra 
were calibrated in wavelength using comparison spectra of 
an HeAr arc taken at the beginning and end of each night, 
and additionally before and after each object. The wave- 
length calibration was done very carefully in order to later 



subtract the sky as accurately as possible. The IRAF rou- 
tine "background" was used to substract the contribution 
of the sky by interpolation of the background detected 
in windows close to and on both sides of the object. Us- 
ing this method, the sky features were successfully sub- 
stractcd. 

2.2. Atmospheric extinction 

The spectra were corrected for atmospheric extinction 
with the aid of mean extinction coefficients for La Silla. 
Molecular absorption bands of O2 (the A band at 7620 A 
and the H2O bands, bands in the region 7100-7450 A and 
8100-8400 A ) were also evident in the spectra. They were 
removed separately from all the observed spectra. These 
absorption bands are composed of many closely spaced ab- 
sorption features, unresolved in our spectra. Many of the 
features which comprise the A band are optically thick 
and therefore nearly independent of zenith distance, but 
most of those in the H2O bands are not independent of 
time and zenith distance. In principle, the best method 
to remove these lines is to observe standard stars as close 
as possible in zenith distance and in time to that of the 
program object in order to obtain the best correction (c.f. 
Osterbrock, Shaw & Vcillcux 1990). But as these authors 
pointed out, the star closest to the objects's observing 
time and zenith distance is not always the one producing 
the best correcting spectrum. 

The procedure we used was the following: we observed 
three different standard stars whose frames were reduced 
in the same way as described before. After modeling the 
atmospheric bands with them, we found that the best re- 
sults were obtained with Feige 110, a dwarf with no in- 
trinsic absorption features in the spectral range observed. 
A 1-dimensional integrated spectrum of the standard star 
was obtained by adding all the light along the spatial di- 
rection. A smooth fit to the continuum was done and then 
the original spectrum was divided by the fit. The result 
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Fig. 1. 

Spectra of observed Radiogalaxies, nuclear (first column) and extended emission line regions (second and third columns). The 
expected position of the Fl line is indicated with | . Flux is given in units of erg s _1 cm~ 2 A -1 . 
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Fig. 2. 

Spectra of observed Liners and Seyfert 2 galaxies. The liner spectra have been extracted from the whole spatial extension of 
the object. The expected position of the Fl line is indicated with | . Flux is given in units of erg s _1 cm~ 2 A -1 . 



was a normalized spectrum which kept the features due to 
atmospheric absorption. Due to the dependence on zenith 
distance and time, we had to model the bands for each 
object, i.e. to construct a specific restoring spectrum for 
each object, multiplying the normalized spectrum by the 
appropriate factors. 

2.3. Flux calibration 

The atmospherically corrected spectra of the standard 
stars were used to obtain the flux calibration. For each 
night, we built a mean response curve from the two stan- 
dard stars observed that night. 

2.4- Template galaxy subtraction 

We checked that effects of stellar absorption features in- 
trinsic to the galaxy could be neglected: the calcium dou- 
blet is on top of the raising part of a molecular band whose 
narrow line components are not resolved. Its only effect is 
to change the slope of the continuum under the doublet 
and we can easily correct for this effect, fitting a smooth 
continuum of appropriate slope in such region. Ha, a line 
which has been used as a reference in our prediction of 
the Fl line flux, can be underestimated due to the under- 
lying Ha absorption line. However, in our objects the Ha 
emission is so strong that the difference is less than 5% 
of the emitted flux. This was estimated by comparing the 
EW of the emission line in our objects with the EW of the 
absorption feature from a template elliptical galaxy. 



2.5. Data analysis 

2.5.1. Extraction of the spectra 

For those objects for which no extended emission lines 
were detected, we compressed the whole spatial extension 
in order to extract a 1-dimcnsional spectrum. For those 
objects showing extended emission lines, we separated it 
into two different components: the nucleus, and the ex- 
tended emission line spectrum. In order to achieve this, 
we subtracted from the spatial profile of the brightest 
spectral line, the spatial profile from the nearby contin- 
uum. In this way, we are left with the spatial profile of 
the line emitting regions. In some cases, emission blobs 
appear well detached from the nuclear emission. In other 
cases, however, the extended and nuclear emission joined 
smoothly. In order to avoid contamination of the EELRs 
by the unresolved nuclear spectrum, we fitted a Gaussian 
(PSF) to the nuclear region. The residuals after subtract- 
ing this Gaussian were considered to pertain to the EELR. 
The resulting spectra are shown in Figs. 1, 2 and 3. 

It is possible that we may be venturing close enough 
to the nucleus to sample higher densities, more typical of 
the classical narrow line region (with densities higher than 
10 3 cm -3 ). To verify this, when it was possible, we used 
the [SII] doublet to estimate the density. The values were 
always below 100 cm~ 3 , which indicates that the spectra 
are dominated by low density gas. In any event, even if 
the NLR dominated the emission, calculations show that 
at high densities the line Fl (in the dustfree case) if any- 
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Fig. 3. 

Spectra ot objects with properties associated with cooling flows. The spectra of PKS1404-267 and A2029 have been extracted 
from the whole spatial extension of the object. The expected position of the Fl line is indicated with | . Flux is given in units 
of erg s _1 cm -2 A -1 . 
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thing becomes tronger than many other lines as a result 
of its very high critical density and of the increase in elec- 
tronic temperature. Our conclusions about inferring dust 
or not are therefore not affected by nuclear contamination. 
(Higher density models would in fact imply a stronger cal- 
cium depletion). 

2.5.2. Line measurements 

IRAF routines were used to measure the emission line 
fluxes. For the blends, decomposition procedures were 
used, fitting several Gaussians at the expected positions 
of the components. The fluxes correspond to the values 
derived from the sum of Gaussians which better fitted the 
profile. 

The upper limits on [Call] emission were measured 
directly from the noise level in the data. 

3. Observed and predicted Fl fluxes 

3.1. Prediction of Fl fluxes 

The second component of the [Call] doublet (F2) is the 
weakest and is blended with the [OII]AA7320,7330 multi- 
plet. The first component is the strongest and lies well 
aside, some 30 A shortward of [Oil]. It is straightfor- 
ward to isolate it given a reasonable spectral resolution. 
Therefore, for simplicity, we have based our study on the 
measurement of this line. The observed values are com- 
pared with the flux that photoionization models predict. 
MAPPINGS (Binette et al. 1993a,b) is the multipurpose 
photoionization-shock code which we used for the model- 
ing. It considers the effects of dust mixed with the ionized 
gas: extinction of the ionizing continuum and of the emis- 
sion lines, scattering by the dust, heating by dust pho- 
toionization and depletion of heavy elements. The appro- 
priate input parameters for the models are justified in Pa- 
per I. 

The predicted value of the Fl flux has been computed 
in the following way: we used as reference line (ref line), 
the strongest and/or the easiest line to measure, like 
Ha, [SII]AA6716,6731, the [OII]AA7320,7330 blend and/or 
[ArIII]A7136. The appropriate photoionization model al- 
lows the prediction of R = rcf F } inc , from which we can 
easily deduct the Fl flux: 

PredictedFl flux = R x Measured flux in the reference 
line. 

In Table 2, we show for each object and for each spa- 
tial region the reference lines used, the measured flux, the 
predicted R = —tt- — ratio, the measured Fl flux and its 

1 rcr. line 7 

theoretical value. For the nuclear regions, we always base 
our predictions on the blend [Oil] AA7320,7330 as refer- 
ence line. The density in the nuclear narrow line region 
spans a wide range up to 10 6 cm~ 3 . The critical density 
of the [Call] forbidden lines is very high, and the flux 
of the lines under the nuclear conditions are predicted to 
be far stronger than the [Oil] lines (Fcrland 1993). The 



low critical density [Oil] line can be de-excited as a result 
of the high nuclear density so that if [Call] remains much 
smaller than [Oil] , it is actually a more stringent test since 
high densities would have actually helped to increase the 
[Call] /[Oil] ratio. To measure the flux in the [Oil] blend, 
we assume that there is a negligible contribution from the 
F2 component. This assumption is reasonable: it is fainter 
that Fl and, in any event, it is near or under the detection 
limit in all cases. We can therefore assume that the [Oil] 
multiplet is not contaminated by F2. 

3.2. Comparison with observations 

The comparison between the measured and predicted Fl 
flux values shown in table 2 demonstrates that, when- 
ever the calculations were possible (sometimes, however, 
there was not any reference line available), it is always 
fainter (non detected in most cases) than expected. This 
result is common to all the regions considered here: nu- 
clear, EELRs and cooling flow filaments. It demonstrates 
that calcium is depleted and, therefore, the gas is mixed 
with dust, both in the nuclear as well as in the extended 
emitting gas in radio galaxies and cooling flow filaments. 

Donahue & Voit applied this same test to the emission 
line nebulae in cluster cooling flows (1993). They inter- 
preted the observed lack of [Call] emission as depletion 
within the ionized filaments of calcium onto dust grains. 
Although for quite a different sample of objects, our work 
support these results, namely that the ionized gas which 
emits the low excitation lines observed in nuclear regions 
of AGNs, in large scale ionized gas in RGs, and in the 
filaments in objects classified as 'cooling flows' must all 
contain dust. 

The most important result is the confirmation that 
dust exists mixed with the ISM in low z radio galaxies. 

4. Discussion 

4-1. Implications on the origin of the gas 

As explained in Paper I, there are some pieces of evidence 
indicative of dust in the extended ionized gas in active 
galaxies, mainly derived from polarization measurements, 
which show the existence of scattered nuclear light over 
large spatial scales, although it is difficult to discriminate 
between dust and electrons as the scattering agent. Our 
results confirm that dust exists mixed with the ISM in low 
redshift radio galaxies. 

As mentioned before, this produces discrepancies with 
the traditional cooling flow theory, which would now be 
required to explain the formation of dust in a shorter 
time than the cooling time! According to this theory (see 
Fabian 1994 for a review), the galaxies, groups and clus- 
ters of galaxies were formed out of gas that collapsed grav- 
itationally. During this process, gravitational energy was 
released which heated the clouds and a hot X-ray atmo- 
sphere resulted. Inhomogeneities in the gas cause matter 
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Table 2. Emission line fluxes for the observed sample of galaxies (in units of erg s~ x cm~ 2 A -1 ). The superscripts M and p 
indicate measured and predicted values respectively. The last two columns show the predicted and measured values for the Fl 
line. For those objects where the calculations were possible, the observed line is clearly fainter than predicted by photoionization 
models. 



AT 

JName 


Spatial reg. 


rej line 


r lux 


Fl V 

ref line 


Flux(Fl) p 


r lux\b 1) 


NGC1052 


Total 


None 








< 6.99e-16 


NGC6215 


Total 


None 








< 3.34e-16 


NGC7552 


Total 


None 








< 3.95e-16 


NGC7714 


Total 


[OILU7325 


1.59e-14 


>1.00 


> 1.59e-14 


< 6.84e-17 


PKS1404-267 


Nucleus 


None 


— 


— 


— 


< 1.79e-16 




5" South 


None 


— 


— 


— 


< 5.23e-17 


PKS2014-55 


Nucleus 


[OII]A7325 


1.54e-15 


>1.00 


>1.54e-15 


< 3.34e-17 




EELR (E) 


Ha 


1.47e-15 


~ 0.45 


6.62e-16 


< 2.60e-17 




EELR (W) 


Ha 


1.06e-15 


~ 0.45 


4.77e-16 


< 3.18e-17 


PKS2152-69 


Nucleus 


[OII]A7325 


7.48e-15 


>1.00 


>7.48e-15 


1.92e-16 




EELR 


[SII]A6731 


1.33e-15 


~ 0.40 


5.32e-16 


< 3.03e-17 




Cloud 


[ArIII]A7136 


3.86e-16 


- 4.50 


1.74e-15 


< 2.24e-17 


PKS2158-380 


Nucleus 


[OII]A7325 


4.03e-15 


>1.00 


>4.03e-15 


3.15e-16 




EELR (E) 


[ArIII]A7136 


7.28e-16 


~ 5.0 


3.64e-15 


< 2.55e-17 




EELR (W) 


[ArIII]A7136 


3.00e-16 


~ 5.0 


1.5e-15 


< 5.00e-17 


PKS2356-61 


Nucleus 


[OII]A7325 


9.39e-16 


>1.00 


>9.39e-16 


2.78e-16 




EELR 


Ha 


5.36e-16 


~ 0.18 


9.6e-16 


< 1.03e-17 


PKS2300-18 


Nucleus 


[OII]A7325 


1.80e-15 


>1.00 


>1.80e-15 


4.37e-16 




EELR 


Ha 


3.61e-16 


~ 0.40 


1.44e-16 


< 1.97e-17 


2A 0335+096 


cD galaxy 


[SII]A6725 


7.46e-16 


~ 0.40 


2.98e-16 


< 2.72e-17 




EELR (E) 


[SII]A6725 


2.24e-15 


~ 0.40 


8.96e-16 


< 1.29e-17 




EELR (W) 


[SII]A6725 


4.22e-16 


~ 0.40 


1.68e-16 


< 1.40e-17 




companion 


[SII]A6725 


7.89e-15 


~ 0.40 


3.16e-15 


< 1.69e-17 


A2029 


Total 


None 








< 1.43e-16 


A2597 


Nucleus 


[OII]A7325 


1.79e-15 


>1.00 


>1.79e-15 


5.14e-16 




EELR 


Ha 


2.32e-15 


~ 0.50 


1.16e-15 


< 5.86e-17 



from the hot atmosphere to cool down and fall towards 
the center of the galaxy or cluster. In this picture, the 
resulting filaments would eventually become visible at op- 
tical wavelengths and would emit strong lines. For massive 
galaxies and cluster, the cooling process would have been 
slower than for normal galaxies and the hot atmosphere 
would still exist, with typical temperatures of several mil- 
lion K, and be emitting strongly in the X-ray band. 

If this is actually the origin of the gas in the EELRs 
(and cooling flow filaments), it should clearly be devoid of 
dust. Any dust introduced in the hot intracluster medium 
would be sputtered and rapidly destroyed (time scale of 
the order of 10 7 years) (Draine & Salpeter 1979), much 
before the filaments cooled down and became visible. Is 
there a way to introduce dust in the filaments during the 
cooling process? The common place where dust is formed 
is stellar atmospheres. Do stars exist in the accreted gas? 
Such gas must be deposited in some form and it has been 
generally thought that the more plausible fate for most 
of the gas is the formation of new stars (e.g., O'Connell 
& McNamara 1989; Fabian, Nulscn & Canizares 1991). 



Although there are some indications of star formation in 
cooling flow galaxies, the latter seems to be taking place in 
the inner parts (over the central few kpc) (Cardiel, Gorgas 
& Aragon-Salamanca 1995; Fabian 1994). There is a lack 
of evidence of star formation (or any kind of stars) at large 
distances, out of the main body of the galaxies, where the 
EELRs filaments are still visible. 

In high z RGs, there is an extended blue continuum 
aligned with the radio axis which has often been inter- 
preted as due to young stars. However, as I mentioned 
already before, further results (e.g., Tadhunter, Fosbury 
& di Serego 1988; Cimatti et al. 1993) proved that this 
continuum is polarized and that the contribution of a hid- 
den quasar continuum scattered by dust and/or electrons 
in the ISM constitutes very likely an important component 
(in fact, this extended polarized UV continuum has been 
used as a proof of the existence of dust in the EELRs of 
high z radio galaxies). This effect has also been detected in 
some radio galaxies at low z (Cimatti & di Serego Alighicri 
S. 1995). Altough young stars may exist, we don't have a 
clear idea about their overall contribution. 
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There is another mechanism which can form dust: if 
most of the cooled gas from a flow does not form stars 
with normal IMF, maybe it remains as cold clouds or as 
low-mass stars. There are evidences of X-ray absorption 
(White et al. 1991; Mushotzky 1992; Allen et al. 1993) 
in cooling flows in cluster of galaxies. The absorbing ma- 
terial could be in the form of cold gas embedded in the 
hot intracluster medium of the cooling flow (White et al. 
1991; Daines, Fabian & Thomas 1994). This cold gas, very 
slightly photoionized by X-rays from the surrounding hot 
corona, can become molecular (Ferland, Fabian & John- 
stone, 1994). Fabian, Nulsen & Canizares (1994) propose 
that the conditions suitable for dust to form may occur in 
this cool gas, through the condensation of gaseous parti- 
cles. The authors propose that the gas which is cooling to- 
wards the center could be a mixture of both the molecular 
gas and the very hot gas and, therefore, could contain the 
existing dust. However, even if such a scheme was possible, 
the dust grains will not remove the calcium which already 
existed in the hot gas; the temperature is too high for the 
condensation of calcium onto the dust grains. Therefore, 
we should observe the Call lines from the hot gas when it 
cools down, even if it contains dust. 

On the other hand, the interaction scenario (gas ac- 
creted from outside the galaxy, as a result of recent tidal 
interactions or mergers between two or more components) 
predicts the existence of dust mixed with the gas, the 
one already existing in the interacting objects. Some mor- 
phological evidences and theoretical ideas (see Paper I) 
support this scene. Heckman et al. (1986) showed that 
a large fraction of powerful radio galaxies have morpho- 
logical features (shells, tails, loops, etc) similar to those 
produced in numerical simulations of galaxy interactions 
(e.g., Toomre and Toomre 1972; Quinn 1984). Kinematic 
measurements show that the radio galaxy EELR generally 
have a high specific angular momentum which is difficult 
to reconcile with the cooling flow picture (Tadhunter, Fos- 
bury & Quinn 1989). 

Therefore, if there is a common origin for the EELRs 
of all radio galaxies, our results suggest that it is mergers 
or tidal interactions. However, we don't exclude the pos- 
sibility of an origin which is not universal, that is, which 
may differ from object to object. 



R.A.E. Fosbury, J. van Loon, R. Pelletier and I. Perez-Fournon 
for useful discussions. 
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5. Conclusions 

We have confirmed that the gas in extended emission line 
regions in radio galaxies at low z is mixed with dust. 

Our results support the existence of dust mixed with 
the gas in the Narrow Line Region and in the cooling flow 
filaments. 

If there is an universal origin for the EELRs, the exis- 
tence of internal dust favours mergers or tidal interactions 
as the most plausible scenario. 
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